1 Detection of Extrasolar Planets by Gravitational 
Microlensing 



David P. Bennett 
University of Notre Dame 

No Institute Given 

Summary. Gravitational microlensing provides a unique window on the properties and 
prevalence of extrasolar planetary systems because of its ability to find low-mass planets 
at separations of a few AU. The early evidence from microlensing indicates that the most 
common type of exoplanet yet detected are the so-called super-Earth planets of ~ 10 
Earth-masses at a separation of a few AU from their host stars. The detection of two such 
planets indicates that roughly one third of stars have such planets in the separation range 
1.5-4 AU, which is about an order of magnitude larger than the prevalence of gas-giant 
planets at these separations. We review the basic physics of the microlensing method, and 
show why this method allows the detection of Earth-mass planets at separations of 2-3 
AU with ground-based observations. We explore the conditions that allow the detection of 
the planetary host stars and allow measurement of planetary orbital parameters. Finally 
we show that a low-cost, space-based microlensing survey can provide a comprehensive 
statistical census of extrasolar planetary systems with sensitivity down to 0.1 Earth-masses 
at separations ranging from 0.5 AU to infinity. 

1.1 Introduction 

The gravitational microlensing method relies upon chances alignments between 
background source stars and foreground stars, which may host planet systems. 
These background source stars serve as sources of light that are used to probe 
the gravitational field of the foreground stars and any planets that they might host. 
The relative motion of the source star and lens system allows the light rays from the 
source to sample different paths through the gravitational field of the foreground 
system, and it is changing total gravitational lens magnification of the source star 
with time that provides the observable gravitational microlensing signal. 

The microlensing method is unique among exoplanet detection methods in a 
number of respects: 

1. The amplitude of planetary microlensing signals is large (typically 10%) and 
is approximately independent of the planetary mass. Instead, the source-lens 
alignment necessary to give a detectable planetary signal depends on the planet- 
star mass ratio, q, and so the probability of a detectable planetary signal scales 
as <~ q. 

2. This scaling of the probability of planet detection with the mass ratio, q, is 
shallower than the sensitivity curves for other methods, so microlensing is more 
sensitive to low-mass planets than other methods that are sensitive to planetary 
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mass. The sensitivity of the microlensing planet search method extends down 
to 0.1M ffi . 

3. Microlensing is most sensitive to planets at orbital separations of 1.5-4AU, 
which corresponds to the vicinity of the Einstein ring radius. This range of sep- 
arations also corresponds to the "snow line" where planet formation is most 
efficient according to the leading core accretion model of planet formation. 
Thus, microlensing complements the Doppler radial velocity and transit meth- 
ods, which are most sensitive to planets in very short period orbits. 

4. Microlensing is the only planet detection method that is sensitive to old, free- 
floating planets, which have been ejected from the gravitational potential well 
of their parent stars through planet-planet scattering. Theory predicts that 
such planets may be quite common, and ground-based microlensing can detect 
free-floating gas giant planets, while a space-based survey is needed to detect 
free-floating terrestrial planets. 

5. Since the microlensing method doesn't rely upon light from the host star in 
order to detect its planets, it can detect planets orbiting unseen stars. This can 
make it difficult to determine the properties of the host stars, but space-based 
follow-up observations can detect the host stars for most planets discovered by 
microlensing. 

6. A space-based microlensing survey would provide a nearly complete statisti- 
cal census of extrasolar planets with masses down to 0.1M® at all separations 
> 0.5 AU. This includes analogs of all the Solar System's planets, except for 
Mercury. 

Gravitational microlensing differs from other extrasolar planet search techniques 
in a number of aspects. It is a purely gravitational method that doesn't rely upon 
on detecting photons from either the planet or its host star 

While most of the known exoplanets have been discovered with the Doppler 
radial velocity method, the early results from the microlensing method indicate that 
cool, super-Earth or sub-Neptune mass planets are more representative of typical 
extrasolar planets than any of the 200+ exoplanets discovered by radial velocities. 



1.2 Gravitational Microlensing Theory 
1.2.1 The Single Lens Case 

The basic physics of gravitational lcnsing depends only a single input from General 
Relativity, the deflection angle, a, for a light ray passing a mass, M, with an impact 
parameter, r: 

4GM 

a = . (1.1) 

c z r 

With the lens geometry shown in Fig. |1.1[ we have 
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in the small angle approximation. If the lens and source are perfectly aligned, the 
two images merge to form a ring of radius 



known as the Einstein ring radius. {Oe is the angular Einstein radius.) We can now 
rewrite the single lens equation as 

n = ^r, (1-4) 

and it has two solutions: r+,_ = 0.5(5 ± ^Jb 2 + 4R|J). The lensed images are also 
magnified, and the magnification of a source of infinitesimal size can be computed 
using area elements obtained by differentiating eq. |1.4| This yields 



A^-l(-$g.±l) , (1.5, 
2 \u\/u 2 + 4 / 

where u = bj Re is the dimcnsionlcss lens-source separation. The total magnification 
of both images is given by 

u 2 + 2 

A = A + + A^ = ^ r ±£=. (1.6) 

For a lens of M = 1M Q , that is half-way to a source in the Galactic center (at 
Ds = 8kpc), we find 



^ = 4.04AU^^4,(1-.), (1.7) 
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where x = Dl/Ds, so Re is similar to the orbital radius of planets in our own 
Solar System. This also implies that 9g ~ lmas. Since the image separation is of 
order ~ 9e, this implies that images will not generally be resolved with virtually all 



planned and future astronomical instruments (with a few exceptions (Delplancke 



et al. 2001 )). On the other hand, if we assume a typical Galactic velocity of v± 
100 km/sec for the relative velocity between the lens star and the line-of-sight to 
the source, then the typical Einstein radius crossing time for a lens in the Galactic 
disk and a bulge source is = Re/v± ~ 2months. Thus, the main observational 
effect for lensing by stars within the Milky Way is the time varying magnification 
instead of the image separation, and this is why it is referred to as microlensing 
instead of lensing. 

The microlensing light curve is generally described by eq. |1.6| with the lens-source 
separation given by 

2 

+ ul , (1.8) 



't-h 



V V *B 

assuming that the relative motion between the lens and the observer-source line-of- 
sight. Thus, a single-lens microlensing light curve is described by three parameters, 
the time of peak magnification, to, the Einstein radius crossing time or width, Ie, 
and the minimum separation, u , which determines the peak magnification. ii is 
the only parameter that affects the intrinsic light curve shape, as shown in Fig. |1.2| 
but only tE constraints the physically interesting parameters of the event: the lens 
mass, M, the lens distance, Dl, and the relative velocity, v±. 




Fig. 1.2. Example microlensing light curves for a point source and a single lens that moves 
with a constant lens velocity with respect to the observer-source line-of-sight. Light curves 
with u = 0.05,0.1,0.2,0.4, and 0.8 are shown. 



The first microlensing events were discovered in 1993 towards the Large Magel- 



lanic Cloud (LMC) by the MACHO Project (Alcock et al. 19931 and towards the 



Galactic bulge by the OGLE Collaboration. (Udalski et al. 1993 1. The early empha- 



sis of microlensing surveys was the search for dark matter in the Milky Way's halo 
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(Paczynski 19861, but this issue has been largely resolved with the demonstration 



that the excess microlensing seen toward the LMC by the MACHO group) Alcock 



et al. 2000b Bennett 2005) requires at most 20% of the Milky Way's dark matter 



in the form of stellar mass objects, while the results of the EROS group (Tisserand 



et al. 2007) suggest that much of this microlensing excess may be caused by stars 



associated with the LMC itself (|SahuJ [1994| , perhaps in the LMC halo ( |Wu[|1994 l. 

With the dark matter microlensing question mostly resolved, the prime focus 
of microlensing observations has shifted to the detection of extrasolar planets. Mi- 
crolensing was first suggested as a method to find planets by Liebes (19641, but as 



Mao & Paczynski ( 1991 1 pointed out, this requires a consideration of multiple lens 
systems. 



1.2.2 Multiple Lens Systems 



The lens equation for a multiple lens system is a straight forward generalization of 
the single lens equation, eq. |1.4[ but with more than one lens mass, we can no longer 
assume that the source, observer and lens system all lie in a single plane. However, 
as long as the distances to the lens and source (Dl and D$) are much larger than 
the extent of the lens system, we can assume that the lens system resides a single 
distance, and define the Einstein radius of the total lens system mass using eq. |1.3| 
So, as before, we will rcscale all the length variables with Re- 

Because we can no longer define a source-lens-observer plane, we must now define 
the lens and source positions in the 2-dimcnsional "lens-plane" perpendicular to the 
line-of-sight and projected to the distance of the lens (or equivalently, we can just 
use angular variables for the positions of the source and lenses on the sky). The 
double-lens system was first solved using two real coordinates for the lens plane 



(Schneider & Weiss 1986[), but the algebra is much simpler if we describe the lens 



plane with complex coordinates following Bourassa et al. (19731; Witt ( 1990J) and 
Rhie ( 1997 ) . The generalization of eq. |1.4| is 



W 



E 



(1.9) 



where w and z are the complex positions of the source and image, respectively, 
and Xi are the complex positions of the lens masses. The individual lens masses are 
represented by e^, which is the mass fraction of the ith lens mass, so that ^ 4 e; = 1. 
The appearance of the complex conjugates in the denominator in the sum on the 



right side of eq. 1.9 is simply a reflection of the fact that the lens deflection is the 
in the direction from the source to the lens with a magnitude of the inverse of that 
distance. With real coordinates, we would express this as the vector difference of 
the positions divided by this vector squared, but with complex coordinates, we can 
divide through by this vector leaving only its complex conjugate in the denominator. 

If we knew the position of the images, z, in eq. 1.9 then it would be trivial to 
solve for the position of the source. But this is the inverse of the problem that we 
will usually want to solve, which is to find the positions of the images based on 
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a known position for the source. However, the solution of this "inverse" problem 



is the basis of the brute-force, ray-shooting method (Schneider & Weiss 19871 for 
solving eq. |1.9| This method involves taking a large grid of points in the "image 



plane" and propagating them back to the source plane using eq. |1.9| This method 
has the advantage that it can handle very complicated lens mass distributions, but 
it is usually not the method of choice for the analysis of microlensing events. 

The most successful method for calculating multi-lens microlensing light curves 
(Bennett & Rhie 19961 involves solving eq. 1.9 for the positions of the point-source 



images and invoking the ray-shooting method only in the vicinity of images that are 
affected by finite-source size effects. For the majority of the light curve, the finite- 
source calculations are not needed, and we can use the point source magnification 
formula. This formula can be derived from the Jacobian determinant of the lens 
equation (and its complex conjugate): 



dw dw dw dw ^ 

dz dz dz dz 



dw 



dz 



where 



dw 
dz 



v. 

^ (z-Xt) 2 



(1.10) 



(1.11) 



Because eq. |1.10| gives the Jacobian determinant of the inverse mapping from the 
image plane to the source plane, the magnification of each image is given by 



A = 



1 



(1.12) 



evaluated at the position of each image. 

The solution of the lens equation, |1.9| is non-trivial. For the case of two lens 
masses, this equation can be embedded into a fifth order polynomial equation in z, 



which can be solved numerically. This equation has either 3 or 5 solutions (Witt 



1990 Rhie 1997) that correspond to solutions of 1.9 which means that a double 



lens system must have either 3 or 5 images depending on the configuration of the 
lens system and the location of the source. For the triple lens case (which is relevant 
for at least one planetary microlensing event), the lens equation can be embedded 
in a rather complicated tenth order polynomial that has 4, 6, 8, or 10 solutions that 



correspond to physical images (Rhie 2002). This tenth order polynomial equation 



can be solved numerically, although it may require extended precision numerical 
calculations in order to avoid serious round-off errors (Bennett et al., in prepara- 
tion). The case of 4 lens masses, has also been investigated (Rhie 20011, but the 
lens equation has not been converted to a polynomial. 

The most important feature of lensing by multiple masses occurs at the locations 
where J = 0. From eq. |1.12| this implies infinite magnification for a point source. 
(The magnification is always finite for the realistic case of a source of finite angular 
size.) For a single lens, J = only occurs at a single point in the source plane, 
the location of the lens mass, but for lens systems with more than one mass, there 
are a set of one or more closed curves with J = 0, known as critical curves. The 
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source positions corresponding to the critical curves are obtained by applying the 
lens equation, |1.9| and they are referred to as caustic curves. When the source 
passes to the interior of a caustic curve, two new images are created, and it is these 
new images that have infinite magnification for the (unphysical) case of a point 
source. The shape of the light curve of a (point) source crossing a caustic has a 
characteristic form: 

F c 0(x — x c ) , ,„ „ „. 

A= / J^ + A nc , (1.13) 

V x ~ x c 

where F c gives the amplitude of the caustic, A nc gives the magnification of the 
images that are not associated with the caustic and 0(x) — 1 for x > and 
0{x) = otherwise, x is the distance perpendicular to the direction of the caustic 
curve, and x c is the location of the caustic curve. Eq. 1.13 is a good approximation 
to the magnification for a point source when the curvature of the caustic curve 
can be neglected. Note, that the singularity in eq. |1.13| is weak enough so that the 
integral of this formula will yield a finite magnification for a finite size source star. 

Caustic crossings that follow the form of eq. |1.13| are often referred to as fold 
caustic crossings, and they have the feature that there is essentially no warning 
that the caustic crossing is imminent when the caustic curve is approached from 
the outside (i.e. x < x c ). This is because the magnification pattern for a fold 
caustic extends only to the interior of the caustic since it involves the magnification 
of images that only exist inside the caustic curve. However, each caustic curve 
also has at least three sharp pointy features, known as cusps, and the magnification 
pattern extends outward from the cusps on a caustic curve. The magnification scales 
as the inverse of the distance to the cusp, just as in the single lens case, eq. |1.6| 

The path of the source with respect to the caustic curves provides the basic 
characteristics of a multiple lens microlensing light curve. Multiple lens light curves 
frequently have features which match the expected A ~ 0(x)x~ 1 ^ 2 shape of a 
caustic crossing or the A ~ r _1 shape of a cusp approach. But, there are additional 
complications, as the strength of a caustic crossing (F c in eq. |1 .13 1 can vary and 
the angular size of the source star can sometimes be larger than the entire caustic 
curve for a planetary microlensing event. 



1.3 Planetary Microlensing Events 

Planetary microlensing events arc a subset of multiple lens events where the mass 
ratio is quite small. Planetary events have light curves that appear quite similar to 
the single lens light curves shown in Fig. |1.2[ but for a brief period of time they 
deviate from the single lens form and display the characteristics of a binary lens 
light curve. We will define a requirement on the mass ratio q = £2/^1 < 0.03 to 



separate planetary microlensing events from stellar binary events following Bond 



et al. (20041, because q ~ 0.03 is the approximate location of the "brown dwarf 
desert" that appears to separate stellar from planetary secondaries. We will also 
initially only consider events with only one detectable planet, as these represent the 
majority of planetary microlensing events and this will simplify the discussion. 



8 Exoplanets via Microlensing 



The caustic structure of a binary lens system is determined by the mass ratio 



and the separation of the lenses (Schneider & Weiss 19861. For a separation d <C 1 



(in units of Re), there are three caustics, two triangular caustics with 3 cusps each 
and a caustic close to the center of mass which has 4 cusps. These merge into a 
single caustic with 5 cusps at d ~ 1, which splits into two caustic curves with 4 
cusps each for d ^> 1. For small values of the mass ratio, q, the division between 
these regimes occurs near d ~ 1, so most events have multiple caustic curves. As 
shown in Fig. |1.3| 



OB05-071 : d=0.759, q=0.006 caustic topology 



OB05-071 : d=1 .299, q=0.006 caustic topology 



Fig. 1.3. Grey-scale magnification patterns for two models of the planetary microlensing 
event OGLE-2005-BLG-71. The darkness of the image is proportional to the logarithm of 
the magnification. The nearly vertical lines in each panel indicate the path of the source 
star for each model. The source passes close to the central, or stellar, caustic, and as 
discussed in § |1.4[ the magnification pattern in the vicinity of the central caustic is similar 
for planetary systems related by d <-> l/d, The magnification patterns for the planetary 
caustics (on the outside of each panel) are clearly very different and easy to distinguish. 
(These images are provide courtesy of Daniel Kubas.) 



1.3.1 Planetary Caustic Perturbations 

There are two classes of caustic curves in planetary microlensing events: planetary 
caustics and the stellar or central caustic. The planetary caustics result when one 
of the two light rays in Fig. |1.1| passes close to the planet and is deflected by the 
planetary gravitational field of the planet. The locations of the caustics are given 
by 

s c = d-l/d. (1.14) 

If we see a planetary deviation at a point where the best fit single lens light curve 
predicts a magnification, A, we can find the corresponding u value by inverting 
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eq. |1.6| to get u c — s c , and then solve for d by inverting eq. |1.14| to give 



d = 



(1.15) 



The two solutions to eq. |1.15| for a given u c value are referred to as major and minor 
image perturbations. They were first studied in detail by Gould & Loeb ( 1992 1 



who showed that many features of the planetary caustics and their magnification 



pattern could be explained by the simpler Chang- Refsdal lens system (Chang & 



Refsdal 1979 1984 1. Several important features of the planetary caustic light curve 



perturbations can be seen in Fig. |1.3| The left side of the left panel of this figure 
shows the magnification pattern of the two roughly triangular minor image caustic, 
which arc generated by a planet with d < 1. As with all caustics, there is excess 
magnification in the interior of the caustic curve, as well as extending outward from 
the cusps. But there is also a very pronounced magnification deficit in between 
the two minor image caustics, where the magnification is substantially below the 
single lens magnification. In contrast, the magnification pattern for the major image 



caustic (shown on the right side of the right panel in Fig. 1.3 is predominantly 



positive, with only small magnification deficits very close to the caustic curve, away 
from the 4 cusps. 



1.3.2 Stellar Caustic Perturbations 

It was originally suspected that the planetary caustic perturbations would be the 
best way to detect planetary signals in microlensing events, but |Grie st & Sa fizadeh| 



( 1998 1 argued that there were a number of advantages to searching for planetary 
light curve perturbations due to the stellar caustic. They showed that the planet 
detection efficiency for each high magnification event was substantially higher than 
for events of more modest magnification. While the higher planet detection effi- 



ciency for higher magnification events was seen in previous work ( Bolatto & Falco 



1994 Bennett & Rhie 19961, Griest & Safizadeh (19981 emphasized that this effect 



is quite dramatic and that this fact could be used to increase the observational 
planet detection efficiency. In the same year that the Griest & Safizadeh paper was 
published, the MPS and MOA Collaborations demonstrated this method with ob- 
servations of the MACHO-98-BLG-35 event. The subsequent analysis showed ( |Rhie| 



et al. 2000 ) that the lens star for this event did not have any Jupiter-mass planets 
with a projected separation of 0.6-8 AU. 

The high planet detection efficiency for high magnification events is particularly 
useful when a large number of microlensing events are discovered by the microlensing 
survey groups. This is the current situation, as the OGLE-3 and MOA-2 surveys 
combine to detect > 700 microlensing events in progress toward the central regions 
of the Milky Way between February and October of each year. Relatively sparse 
monitoring of events (i.e. one or two observations per day) is required to predict 
most high magnification events in advance, and this allows observing resources can 
be focused on events with a high planet detection efficiency. 
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One important consequence of the high planet detection efficiency for high mag- 
nification events is that the chances of detecting multiple planets in such events 



are greatly enhanced (Gaudi et al. 19981. Indeed, the first multi-planet system dis- 
covered by microlensing is shown below in § |1.5.2| There is, however, a potential 
downside to this higher sensitivity to multiple lens masses. The signals for all the 
detectable lens masses will be concentrated in the very high magnification part of 
the light curve, and this could make it difficult to work out the details of multiple 
planet systems that are detected in microlensing events. Thus, the development of 
efficient light curve modeling methods for lens systems with three or more masses 
is an important active area of current research. 

A final advantage of high magnification events is that they allow planet detection 
with relatively faint source stars. This makes it much easier to detect the planetary 



host star with follow-up observations (Bennett et al. 2007a) as explained in 



1.4.1 



1.3.3 Finite Source Effects 

Microlensing is arguably the exoplanet search technique that is most sensitive to 
low-mass planets, and lower limit in sensitivity of the microlensing method is set 
by the finite angular size of the source stars. Roughly speaking, when the angular 
radius of the planetary Einstein ring, y/ep0Ei is much smaller than the source star 
angular radius, 8*, we expect that planetary signal to be washed out. But this is 
only a crude, order-of-magnitude estimate, and a full finite source solution to the 
lens equation, |1.9| is required to determine the precise limits on the microlensing 
planet detection method set by the finite angular size of the sources. 

Full finite source planetary microlensing light curves were first calculated by 



Bennett & Rhie (19961 using the methods described in § 1.2.2 Results of these 
calculations are reproduced in Figs. |1.4| and |1.5| Fig. |1.4| shows a series of planetary 
light curves with planetary mass fractions of e — 10~ 4 and 10~ 5 . For a typical lens 
star mass of ~ 0.3Af Q , these correspond to lAf® and 10M®, respectively. The finite 
source light curves are characterized by the source star radius in Einstein ring units: 



0*/9e- The p values shown in Fig. 1.4 are 0.003, 0.006, 0.013, and 0.03, and 



these span the expected range of p for a low mass planetary host star in the Galactic 
bulge with a source star ranging in radius from IRq to 10-Rq, which is a typical 
radius for a "red clump" K-giant in the bulge. A number of the planet detections to 
date actually have p values in the 0.4-1 x 10~ 3 range because the lens stars reside in 
the bulge and have a larger than average mass. (Both of these imply a larger 8e-) 
Several general trends are apparent from Fig. |1.4| First, the planetary deviations 
are easily detectable for p = 0.003, but the signals are much weaker for p = 0.03. 
This implies that lAf® planets are easily detected with main sequence source stars, 
but planets of 10M® are close to the lower limit of detectability for giant source 
stars. Another notable feature of these light curves is that the planetary signals 
with d = 0.8 are more easily washed out by the finite source effects than those with 
d = 1.3. This is a consequence of the large magnification deficit seen between the two 
planetary caustics for a minor image perturbation, as shown in the left hand panel of 
Fig. 1 1 - 3 1 When a large finite source effectively averages over the vicinity of a minor 
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Fig. 1.4. Microlensing lightcurves which show planetary deviations are plotted for a mass 
ratio of e = 10 -4 & 10 -5 and separations of d = 1.3 & 0.8. The main plots are for a stellar 
radius of p — 0.003 while the insets show light curves for radii of 0.006, 0.013, and 0.03 
as well. The dashed curves are the unperturbed single lens lightcurves. For each of these 
lightcurves, the source trajectory is at an angle of sin -1 0.6 with respect to the star-planet 
axis. The impact parameter uo = 0.27 for the d — 0.8 plots and uo = 0.32 for the d — 1.3 
plots. 



image planetary caustic, the positive and negative deviations effectively cancel each 
other out (Bennett & Rhie 1996 Gould & Gaucherel 1997). In contrast, the major 
image planetary caustic magnification deviation pattern is mostly positive, so the 
finite source effect merely smoothes it out. 

Fig. |1.5| shows how the planet detection probability varies as a function of d 
for the same e and p values used for Fig. |1.4| The greater tolerance of deviations 
with d > 1 to finite source effe cts i s apparent for e = 10~ 4 , p = 0.03 and e = 10~ 5 , 
p = 0.013. The behavior of Fig. *~ 



1.5 



near d = 1 is even more interesting. For e = 10 
and p = 0.003 or 0.006, the detection probability reaches a maximum at d ~ 1, but 
for p = 0.013, the probability has a local minimum at d = 1, and for p = 0.03 or 
any of the p values with e = 10~ 5 , the detection probability = at d = 1. This 
is due to the fact that for d « 1 the planetary and stellar caustics merge for form 
a relatively large single caustic that is extended along the lens axis. This caustic 
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Fig. 1.5. The planetary deviation detection probability is plotted for different values of 
the planetary mass ratio, e, and the stellar radii, p (in units of Re- A planet is considered 
to be "detected" if the lightcurve deviates from the standard point lens lightcurve by more 
than 4% for a duration of more than tE /200. 



is large, but relatively weak, and it has associated positive and negative deviation 
regions that tend to cancel when averaged over by a moderately large finite source. 

However, some features of Fig. 1.5 are dependent on the somewhat arbitrary 
choice of the event detection threshold, and the sensitivity of a real observing 
strategy can differ from this. In fact, the planetary deviation detected in event 
OGLE-2005-BLG-169 flGould et ah] |2006[ ) would not have passed the selection cri- 
teria for Fig. 1.5 but the planet is nevertheless detected with a strong signal. The 
reason for this is that it was identified as a very high magnification event with a 
very high sensitivity to planets, and for this reason it was observed much more 
frequently than most events with potential planetary signals. The additional ob- 
servations provided enough additional signal to allow the definitive detection of a 
relatively low-amplitude signal. 

As a practical matter, finite source effects imply a lower planetary mass limit 
of M p £ 5M0 for giant source stars in the bulge, and a limit of M p £ 0.05M® for 
bulge main sequence stars. Thus, searches for terrestrial exoplanets must focus on 
main sequence source stars. 



1.4 Planetary Parameters from Microlensing Events 

The determination of the properties of the lens systems that are detected in mi- 
crolensing events is often a serious challenge. The simple form of the microlensing 
light curves shown in Fig. |1.2| is an advantage when trying to identify microlensing 
events, but as I mentioned in § |1.2.1| in a single lens event, it can also be a drawback 
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when trying to interpret observed microlensing events. For most single lens events, 
it is only the ig parameter that constraints the physically interesting parameters 
of the event: the lens mass, M, the lens distance, Dl, and the relative velocity, v±. 
The single lens parameters u$ and to don't constrain lens system parameters that 
are of much interest . 

In addition to the parameters needed to describe a single lens event, a planetary 
microlensing event must have three additional binary lens parameters: the planetary 
mass ratio, q = e/(l — e), the star-planet separation, d, (which is in units of Re), and 
the angle between the star-planet axis and the trajectory of the source with respect 
to the lens system, 6. So, two of these new parameters, q and d, directly constrain 
planetary parameters of interest, although d is normalized to Re, which may not be 
known. Most planetary light curves, at least those for low-mass planets, also have 
caustic crossings or a close approach to a cusp that reveal light curve features due 
to the finite size of the source star. This enables the source radius crossing time, t*, 
to be measured. 

The determination of the star-planet separation and the planetary mass frac- 
tion is usually quite straightforward from the microlensing light curve. For events 
at moderate magnification, due to the planetary caustic, the separation can be de- 
termined by the magnification predicted by the single lens model that describes 
the event outside the region of the planetary deviation following eq. |1.15| This still 
leaves an ambiguity between the d < 1 and d > 1 solutions, but this is easily re- 
solved by the drastically different magnification patterns in the vicinity of major 
image and minor image caustics, as shown in Fig. |1.3[ The planetary mass fraction, 
q, can generally be determined by the duration of the planetary perturbation. In 
some cases, if the time scale of the deviation is similar to or smaller than t*, both 
q and i* determine the deviation time scale, but good light curve coverage with 



moderately precise photometry allow both q and t* to be determined (Gaudi & 



Gould 19971 



The situation is somewhat different for high magnification, stellar caustic de- 
viation events. Dominik (19991 pointed out an approximate degeneracy in the 
properties of the stellar caustic under the transformation d — > 1/d, which means 
that there may be a d <-*■ 1/d ambiguity in the modeling of stellar caustic planetary 
events. This is apparent from the magnification patterns shown in Fig. |1.3| For as 
source trajectory nearly parallel to the lens axis or for d ~ 1, this degeneracy breaks 
down, so the ambiguity disappears. With precise photometry it is usually possible 
to to distinguish between the d < 1 and d > 1 solutions, and this has been the case 
for all events observed to date. 



1.4.1 Angular Einstein Radius 

A large fraction of planetary light curve deviations exhibit finite source effects that 
allow the source radius crossing time, t* , to be measured. This is the case for most 
detectable events with a planetary mass, M p but for gas giant planets 

of M p ^ 3OOM0, it is possible to detect a planetary deviation without the source 
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crossing a caustic or coming close enough to a cusp to display finite source effects. 
So, is measurable for most, but not all planetary microlensing events. 

When is measured, it is possible to place an additional constraint, as long 
as the angular radius of the source star, 9* , can be estimated, because the angular 
Einstein radius is given by 

e E = ^ . (i.i6) 

The angular radius of the source star can be measured if the brightness and color 
of the source are known with the use of empirical color-angular radius relations 
van Belle 1999 Kervella et al. 20041. In the crowded fields where microlensing 



events are observed, the most reliable measure of the source star brightness and 
color comes from the light curve models, which include the source brightness as a 
model parameter. So, although it is sufficient to measure the detailed light curve 
shape in a single passband, it is important to obtain a few measurements during 
the microlensing event in at least one additional passband so that the light curve 
fit will also reveal the color of the source. It is also important to estimate the 
extinction towards the source. With measurements in only two colors, such as V 
and J, the extinction can be estimated by comparison to the red clump giant stars 
within an arc minute or two of the target star (Yoo et al. 2004). While this does 



not yield a precise measure of the extinction to the source, note that an error in 
the extinction to the source will affect both the estimated intrinsic brightness and 
color of the source. Fortunately, the extinction-induced brightness and color errors 
have the opposite effect on the estimated source star radius. This partial cancelation 
implies that the estimated value is not very sensitive to to the uncertainty in the 
extinction. 

A more precise estimate of 9„ can be obtained with observations during the 
microlensing event in more than two passbands, particularly if one of the passbands 
is in the infrared because the optical-IR color-radius relations are much more precise 



than the optical ones (Kervella et al. 2004) and because extinction is much lower in 



the IR than in the optical. Observations in 3 or more colors also allow an estimate 
of extinction that doesn't depend the nearby clump giants, with the use of empirical 



color-color relations (Bessell & Brett 19981. 



When the angular Einstein radius is measured, we have the following relation, 



M L = 



c 2 2 D s Dl 
4G E D S -D L 



(1.17) 



which can be considered to be a mass-distance relation because D$ is generally 
known (approximately) from the brightness and color of the source. (The high den- 
sity of stars in the Galactic bulge means that the source is almost always a bulge 



star.) Eq. 1.17 provides a one-parameter family of solutions to the microlensing 
event, and this can be converted to a measurement of the planetary host star prop- 
erties with one additional piece of information. Since the brightness of the source 
star can be determined by the light curve fit, the brightness of the lens star can be 
determined with an image that has sufficient angular resolution to resolve the source 
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Fig. 1.6. The predicted fractional brightness, /i cns = F L /(F S + F L ), of the OGLE-2003- 
BLG-169 lens is plotted in the top panel as a function of mass in the BVIJH passbands. 
(Fl and Fs and the detected flux of the lens and source, respectively. The predicted offsets 
of the centroids of the blended source+lens images in different passbands are shown in the 
bottom panel, assuming that the images are taken 2.4 years after peak magnification. 



and lens stars from the unrelated stars in the field. This generally requires space- 
based imaging with the Hubble Space Telescope (HST), or possibly ground-based 
adaptive optics imaging because of the extreme crowding in the Galactic bulge fields 
where microlensing events are most easily found. (The lens-source relative proper 
motion has typical value /i re i ~ 5mas/yr, so the lens and source are not typically 
resolved from each other until a decade or more after the event.) If the combined 
lens-plus-source image is significantly brighter than the brightness of the source 
from the microlensing fit, then the difference determines the brightness of the lens. 
This then allows the mass of the planetary host (lens) star to be determined using 
a main sequence star mass-luminosity relation (Bennett et al. 2007a I. 



The top panel of Fig. |1.6| shows the predicted brightness of the lens for the 
OGLE-2005-BLG-169 event in the BVIJH passbands. This indicates that the lens 
star will easily be detected if it is a main sequence star, since even a 0.08M Q lens 
star will contribute /t 40% of the ff-band flux and ^ 10% of the /-band flux. 
This case is more favorable than most because of a relatively large 6e value, but 
most cases, the lens star will be detectable in the ii-band unless it is a late M-dwarf 
located in the bulge. However, for Galactic disk lenses at a certain range of distances 
(corresponding to 0.2M o < M < 0.4 for OGLE-2005-BLG-169 in the IJH-bands) 
the mass-distance relation, eq. |1.17| combines with the mass-luminosity relation to 
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yield a nearly flat mass-brightness relation for the planetary host star. In these cases, 
it is useful to have images in shorter wavelength bands, such as V and B because 
this cancelation does generally not occur in the optical and infrared passbands for 
the same range of lens star masses. 




Fig. 1.7. The top-left panel shows the fraction of the source+lens flux for event OGLE- 
2003-BLG-235/MOA-2003-BLG-53 that is predicted to come from the lens in the HST-/, 
V, and B passbands as a function of lens mass. The bottom-left panel shows the predicted 
color-dependent centroid shifts as a function of mass for 1.78 years of relative proper 
motion at ^t r ei = 3.3mas/yr. The measured values of /i en s in the I-band and the color 
dependent centroid shifts and error bars are indicated with their error bars. These are 
plotted at an arbitrary value for the stellar mass (M t ). The centroids of the source+lens 
star blended images in the individual HST/ACS/HRC images are shown in the right panel 
as red circles (7), green squares (V), and blue triangles (B). The crossed error bars are 
the average centroid in each passband. 



High resolution images in multiple colors also allow an independent method 
for estimating the lens star brightness, as shown in the bottom panel of Fig. |1.6 
and the bottom-left panel of Fig. |1.7| Because the lens and source stars usually 
have different colors, the centroid of the blended source+lens image will usually 
be color dependent. So, an additional constraint on the lens star is obtained by 
measuring the centroid offset between the centroids of the blended source+lens in 
different passbands. As indicated in Fig. |1.7[ this effect wa s marginally detected for 



the first planet detected by microlensing (Bennett et al. 2006} with HST images 



taken only 1.8 years after peak magnification. Also, because this color dependent 
centroid shift depends on the relative lens-source proper motion, /i ro i, it can be used 
to help determine 6e for planetary events with no finite source effects, and hence, 
no measurement of . 

The stable point-spread function (PSF) of space-based telescopes, such as HST, 
allows the measurement of the image elongation due to the growing separation of 
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Simulated HST images: 




raw image PSF subtracted binned 



Fig. 1.8. Simulated image stacks of multiple dithered exposures of the OGLE-2005- 
BLG-169 source and lens star 2.4 years after peak magnification using the HST/ACS High 
Resolution Camera (HRC) in the F814W filter band. The top row of images assumed 
a host star mass of M, = 0.08Mq, the middle row assumes M, = 0.35Mq, and the 
bottom row assumes M* = 0.63Mq. In each row, the image on the left shows the raw 
image stack sampled at one half the native HRC (28 mas) pixel size. The central column 
shows the residuals after subtraction of the best fit PSF model, showing the blended image 
elongation along the i-axis due to the lens-source separation. The right hand column shows 
these residuals rebinned to the 28 mas native pixel scale. 

the lens and source stars after the microlensing event. Simulations of this effect 
for the OGLE-2005-BLG-169 event are shown in Fig. 1 1 . 8 1 for three different cases: 
M L = O.O8M , M L = 0.35M Q , and M L = O.63M . This event has a higher relative 
proper motion than most events, but this simulation assume images taken only 2.4 
years after peak magnification, so for other events, it may be necessary to obtain 
the follow-up space-based images ~ 4 years after peak magnification. 

When /i le i can be measured from image elongation and/or the color dependent 
ccntroid shift, then the angular Einstein radius can be determined via 

0E = VreltE , (1.18) 

so that mass-distance relation can be determined even when t* cannot be measured. 
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1.4.2 Microlensing Parallax 



Another method to "solve" a microlensing event involves the microlensing parallax 
effect. This refers to measurements of r e = ReDs/(Ds — Dl), the Einstein radius 
projected to the position of the Solar System. f e can be measured with the help 
of observations of microlensing events by observers at different locations. Because 
the Einstein radius is typically of order Re ~ 1 AU, the observers must generally 
be separated by a distance of ~ 1 AU. The conceptually simplest way to do this is 



to observe an event simultaneously with a satellite in a heliocentric orbit, (Rcfsdal 



1966 Gould 19921, as has recently been done with Spitzcr (Dong et al. 2007) 



However, it is much more common to use the orbital motion of the Earth the measure 



the microlensing parallax effect ( Alcock et al. 1995 Mao 1999 Smith et al. 2002 



Bennett et al. 2002 Mao et al. 2002 1, but for events of very high magnification, it 



is possible to measure this effect with observations from different observatories on 



Earth ( Gould 1997 ) , as has recently been done by Gould et al. (2007, in preparation) . 



A potential complication with this method is that the orbital motion of the 
source star can mimic the effect of the orbital motion of the Earth, but if the signal 
is strong, it is generally possible to detect the characteristic features of the Earth's 
orbit (Poindexter et al. 20051. Another potential complication is that for events 



with tE <C lyr, it is often possible to measure only a single component of the 
two-dimensional f e vector (Smith et al. 2003 1. But, for events with detectable lens 
stars, the two-dimensional relative proper motion, /x rol , it is possible to determine 
the full ? e vector because f e || [i Tel . 

When f e and 9e are both measured, the lens system mass is given by 



Mr 



c 2 . 



AG 



r e VE 



(1.19) 



This method has been used to determine the lens mass for a binary star lens sys- 



tem towards the Galactic bulge (An et al. 20021 and a low-mass M-dwarf lens 



towards the Large Magellanic Cloud (Gould, Bennett, & Alves 20041. The first use 



of this method in a planetary microlensing event is the case of the double planet 
event OGLE-2006-BLG-109, to be published later this year (Gaudi et al. 2007, in 
preparation; Bennett et al. 2007, in preparation). 



1.4.3 Planetary Orbits 

The final property of a planetary system that can be measured is the orbital motion 
of the planet with respect to the star. This is a lower order effect than microlensing 
parallax because we see the effects of both the planet and the star in the light curve. 
So, we are sensitive to the relative velocity between the star and planet, whereas 
the velocity of the Earth around the Sun cannot be separated from the lens-source 
relative velocity. However, the time scale of the planetary deviation is generally 
only a small fraction of the microlensing light curve, and this limits the amount of 
time over which we can detect the orbital motion effects. Also, the typical orbital 
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period of a planet detected by microlensing is ~ lOyrs, so the orbital velocities are 
generally lower than that of the Earth around the Sun. 

The for a planetary deviation of duration At and an orbital period, P, the orbital 
motion during the planetary deviation causes a shift in the planetary lens position 
with respect to the source of order 

Au « At 2 ^ w 0.002 - 0.02 , (1.20) 

assuming a planetary deviation duration of 1 — 10 days. In oder to determine whether 
eq. 1.20 indicates that the effect of orbital motion is detectable, we need to know 
what value of Au is measurable. On thing that limits our resolution in Au is the 
finite angular size of the source star. The typical angular size for a bulge main 
sequence source is 6* ~ 0.5 /^as and a typical angular Einstein radius for a bulge 
event is Be ~ 0.5 mas, so the source radius is typically of order p = 9*/0e ~ 0.001. 
So, if we can detect Au as small as O.lp, then we could be sensitive to Au ~ 10~ 4 . 

In practice, it can be difficult to do this well in the measurement of orbital effects 
because changes in other model parameters can often compensate for the change in 
Au due to orbital motion. In order to retain a constraint on the orbital motion, it is 
generally necessary to have a relatively complicated planetary deviation with more 
than a single caustic crossing or cusp passage that is well sampled by the data. 
Finally, for events with relatively long planetary signals, the orbital acceleration 
can be as large as Au w (At2ir / P) 2 w 4 x 10~ 4 . So, with a very well sampled 
planetary deviation it is also possible to measure the orbital acceleration, as well as 
the velocity. 



1.5 Observational Programs 



There are a variety of different observing programs that contribute to the detection 
of planets via gravitational microlensing. The most basic requirement is to be able 
to identify microlensing events, as was first done by the MACHO Collaboration 



towards the LMC ( Alcock et al. 1993 I and OGLE group toward the Galactic bulge 



(Udalski et al. 19931. Because microlensing observing programs do not yet have 
the resources to observe £ 10 square degrees of the Galactic bulge several times 
per hour, it has been necessary to follow a strategy first suggested by |Gould fc] 
Loeb (1992). Stellar microlensing events must be identified in progress, and then 
followed with a global network of telescopes on an ~ hourly time scale. The MACHO 



(Alcock et al. 1994 19961 and OGLE (Udalski et al. 1994) groups developed real 



time microlensing detection systems within a year after the first microlensing events 
were discovered, and this led to the first spectroscopic confirmation of a microlensing 
event (Benetti et al. 19951. The MOA group began real time detections in 2000 



(Bond et al. 2001 ) and was the first group to employ real time event detection with 



the more advanced difference imaging photometry method (Bond et al. , 2002a ) 



The first microlensing follow-up projects were the Probing Lensing Anomalies 



NETwork or PLANET group (Albrow et al. 19981 and the Global Microlensing 
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Alert Network, or GMAN, (Pratt et al. 1995), which both began taking data in 
1995. The PLANET team followed the |Gould fc Loeb| fll992| strategy, but the 
GMAN group focused more on non-planetary microlensing. A second follow-up 
group focused on exoplanets, the Microlensing Planet Search (MPS) collaboration 

" [1999} , but MPS merged with PLANET in 2004. The 



began in 1997 (Rhie et al. 



final microlensing follow-up group is the Microlensing Follow-up Network or Mi- 



croFUN (Yoo et al. 20041, which began observations in 2003. MicroFUN does not 



follow the Gould & Loeb ( 1992 1 strategy, but instead focuses on high magnification 
microlensing events as suggested by Griest & Safizadeh (19981. 



1.5.1 Early Observational Results 



The most definitive of the early planetary microlensing observational results in- 
volved limits on the presence of planets based on the lack of detection of planetary 
signals. The MPS and MOA groups reported the first planetary limits from a high 
magnification event ( |Rhie et ah 20001. This, was the first demonstration of sensi- 



tivity to Earth-mass planets by any method, except for pulsar timing (Wolszczan 



& Frail 19921. The PLANET group followed with limits from a lower magnifica- 



tion event ( Albrow et al. 2000b) and then a systematic analysis of five years worth 
of null detections (Albrow et al. 2001 Gaudi et al. 2002). They found that less 



than 33% of the lens stars in the inner Galactic disk and bulge have companions 
of a Jupiter mass or greater between 1.5 and 4 AU. These papers claim that their 
limits apply to Galactic bulge M-dwarfs, but this summary of the PLANET result 
neglects an important bias in the events that have been searched for planets. The 
microlensing teams are more efficient at finding long time scale microlensing events 



(Alcock et al. 




2000a 




Sumi et al. 


Hamadache et al. 


2006) The long 



2003 Popowski et al. 2005 Sumi et al. 2006 



to peak magnification, so they can be more efficiently searched for planetary signals. 
As a result, the median time scale of the events search for planets in |Gaudi et al.| 
( |2002 1 is [ts) — 37 days, while the actual efficiency corrected median time scale is 
(tE) = 16 days. This implies that the events that have been searched for planets 
have more massive lens stars and are more likely to reside in the disk than the 
typical Galactic bulge microlensing event. Thus, it is probably the case that most 



of the events searched by |Gaudi et al. ( 2002 1 have lens stars that are more massive 
than an M-dwarf, reside in the Galactic disk, or both. 

In addition to these upper limits on the planetary frequency, there were also a 



number of less-than-certain planet detections. |Rhie fc Bennett ( 1996[ ) showed that 
the very first microlensing event discovered showed a light curve feature that could 
be explained by a planet, but there was a near equal mass binary lens fit that could 



also explain the data. The MACHO group ( jBennett et al. 1997) pointed out that 



there is a good chance that event MACHO-95-BLG-3 was caused by a free-floating 
Jupiter-mass planet. The MPS group found that their data for MACHO-97-BLG-41 



was best explained by a Jupiter-mass planet orbiting a binary star system (^Bennett 



et al. 1999), but the PLANET data for this event favored an orbiting binary star 



interpretation (Albrow et al. 2000a). (Some of the MPS data are now known to 
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be contaminated by moonlight reflecting off the telescope optics.) An analysis by 
the MO A group QBond et alT] |2002b[) showed that the combined MACHO, MO A, 



MPS, and PLANET data for MACHO-98-BLG-35 was consistent with the low S/N 



detection of a terrestrial planet. Finally, Jaroszynski & Paczyhski (20021 showed 



that the event OGLE-2002-BLG-55 had a signal consistent with a planet detection, 
but Gaudi & Han (20041 pointed out that there were other possible explanations. 



Table 1.1. Exoplanets Discovered by Microlensing 



Event Name 



Star Mass Planet Mass Semi-Major Axis Lead Group 



0.63 



+0.07 



OGLE-2003-BLG-235Lb/ 
MOA-2003-BLG-53Lb 
OGLE-2005-BLG-71Lb 0.46 ± 0.04M, 

OGLE-2005-BLG-390Lb 

OGLE-2005-BLG-169Lb 

OGLE-2006-BLG-109Lb 
OGLE-2006-BLG-109Lc 



M„ 830+?™M, 



190 



0.22 
0.49 



+0.21 



+0.14 



0.18 '"S 

0.50 ± 0.05M Q 
0.50 ± 0.05M Q 



1100 ± 100M a 

13; 

226 ± 25M® 
86 ± 10M ffi 



; +5.5 
-2.7 J 



4.3+^AU 



4.4 ± 1.8 AU 

2.6+^AU 

3.2+^AU 

2.3 ± 0.2AU 
4.6 ± 0.5AU 



MOA 

OGLE 

PLANET 

MicroFUN 

MicroFUN 
MicroFUN 



1.5.2 Microlensing Planet Detections 



Table |1.1| summarizes the properties of the planets discovered by microlensing to 



date, including four published microlensing exoplanet discoveries (Bond et al. 2004 



Udalski et al. 2005 Beaulieu et al. 2006 Gould et al. 2006 ) plus a 2-planet system 



that will soon be published (Gaudi et al. 2008, Bennett et al. 2009, in preparation). 
The microlensing discoveries are compared to other known exoplanets in Fig. |1.9| 

The first planet discovered by microlensing is shown in Fig. 1.10| The light 
curve is plotted in units of the source star flux, which is determined by the best 
microlensing model to the event, because the star field is too crowded to determine 
the unmagnificd stellar flux directly. This event was first discovered by the OGLE 
group and announced via their "early warning system" as event OGLE-2003-BLG- 
235 on 2003 June 22. On 2003 July 21, the alert system of the MOA-I microlensing 
survey detected this event and reported it as MOA-2003-BLG-53. The MOA de- 
tection came later because the MOA-I telescope had only a 0.61m aperture and 
has worse seeing conditions than are typical at the 1.3 m OGLE telescope in Chile. 
However, the MOA telescope had a larger field-of-view (FOV), and this enabled 
them to image each of their survey fields ~ 5 times per clear night. As a result, 
MOA was able to detect the second caustic crossing for this event, and arrange for 
the additional observations that caught the caustic crossing endpoint (thanks to 
first author, Ian Bond, who was monitoring the photometry in real time). 

The naming convention for planets discovered is that the name from the first 
team to find the microlensing event is used for the event, so in this case OGLE- 
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Exoplanet Discovery Potential 




semi — major axis (AU) 



Fig. 1.9. The sensitivity of various exoplanet detection methods is plotted in the mass vs. 
semi-major axis plane. Doppler radial velocity detections are shown in black, with 1-sided 
error bars for the m sin i uncertainty. Planets first detected by transits are shown in blue, 
and the microlensing planet discoveries are shown in red. The gold, cyan and green shaded 
regions show the sensitivity of the radial velocity method and NASA's Kepler and SIM 
missions, respectively. The light red and red curves show the sensitivity of current and 
future microlensing planet search programs, and the purple curve gives the sensitivity of 
the proposed Microlensing Planet Finder (MPF) mission. 



2003-BLG-235 takes precedence over MOA-2003-BLG-53. When referring to the 
lens system, we add a suffix "L", and when referring to the source, we add an "S". 
For a lens or source system that is multiple, we add an additional capital letter 
suffix for a stellar mass object or a lower case letter for a planetary mass compan- 
ion. So, OGLE-2006-BLG-109LA, OGLE-2006-BLG-109Lb, and OGLE-2006-BLG- 
109Lc, refer to the star and two known planets of the OGLE-2006-BLG-109 lens 
system. This convention provides names for multiple components of the source star 
system. For example, OGLE-2022-BLG-876Sb would refer to a planetary compan- 



ion to the source star (which would be difficult, but not impossible (Graff & Gaudi 



2000 Lewis 2001 )to detect. 
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Fig. 1.10. The OGLE-2003-BLG-235/MOA-2003-BLG-53 light curve with OGLE data 
in red and MOA data in blue. The top-left panel presents the complete data set during 
2003 (main panel) and the 20012003 OGLE data (inset). The median errors in the OGLE 
and MOA points are indicated in the legend. The bottom panel is the same as the top 
panel, but with the MOA data grouped in 1 day bins, except for the caustic crossing 
nights, and with the inset showing MOA photometry during 20002003. The binary- and 
single-lens fits are indicated by the solid black and cyan dashed curves, respectively. The 
right panel shows the light curve and models during caustic traverse. These models are the 
single-lens case (cyan, long-dashed curve), the best binary lens with q ^ 0.03 (magenta, 
short-dashed line), the planetary lens with caustic entry before day 2835 (green, dotted 
line), and the best overall fit with q — 0.0039(black, solid line). The insets show the second 
caustic crossing and a region of the declining part of the light curve where the best-fit 
nonplanetary binary-lens model fails to fit the data. MOA data on days other than the 
caustic entry and exit (days 2835 ± 0.5 and 2842 ± 0.5) are placed in 1 day bins. 



It is interesting to note that this event was discovered by a procedure that differs 



from both the alert-plus- followup strategy suggested by Gould & Loeb ( 1992 ) and 
the high magnification strategy suggested by Griest & Safizadeh (19981. Instead, 



the planetary deviation was detected in the observations of one of the survey teams, 
and identified in time to obtain additional data to confirm the planetary nature of 
the light curve deviation. We will return to this strategy later in the discussion of 
future microlensing projects given in § |1.6| 

Another notable feature of this event is that the lens star has been identified in 
HST images (Bennett et al. 2006 ). As indicated in Fig. |1.7[ there is an extra source 
of light superimposed at the location of the source star. This is very likely to be the 
lens star, and if so, the HST photometry implies that a fraction, /i cns = 0.18 ±0.05, 
of the total source plus lens flux comes from the lens. During the microlensing 
event, the lens and source were separated by < 0.1 mas, but by the time of the HST 
images, At — 1.78 years after peak magnification, the lens-source separation should 
have grown to Atfi lc \ = 5.9 ± 0.7mas. (/z re i = 3.3 ± 0.4mas/yr was determined 
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Fig. 1.11. Bayesian probability densities for the properties of the planet, OGLE-2003- 
BLG-235Lb, and its host star if it is a main sequence star, (a) The masses of the lens star 
and its planet (M* and M p respectively), (b) the separation, (c) their distance from the 
observer (-Dl); and (d) the I-band brightness of the host star. The dashed vertical lines 
indicate the medians, and the shading indicates the central 68.3% and 95.4% confidence 
intervals. All estimates follow from a Bayesian analysis assuming a standard model for the 
disk and bulge population of the Milky Way, the stellar mass function of Bennett & Rhie 
(2002). 



from eq. 1.16 with input parameters from the light curve model.) This separation, 
plus the mass-distance relation, eq. 1.17| enable to derivation of the curves shown 
in the bottom left panel of Fig. |1.7[ These show the amplitude for the offset of the 
centroids of the blended lens plus source images in different color bands. The HST 
data indicate a marginal detection of this color-dependent centroid shift at a level 
consistent with the assumption that the excess flux is due to the lens. 

With this marginal detection of the color-dependent centroid shift, we can't be 
absolutely sure that the lens star has been detected because it is possible that the 
excess flux could be due to a companion to the source star. It is straight forward to 



deal with this uncertainty with a Bayesian analysis (Bennett et al. 2006 ), and the re- 



sults of such an analysis are shown in Fig. |1.11| The resulting most likely parameter 



values for the event parameters are a host star mass of M* = O.63^q'q 9 M , a planet 
mass of M p = 2.6*g'gMj up , and an orbital semi-major axis of a = 4.3 AU. The 



-0.6 

distance to the lens system is Dl 

j i = 2i.4±g : «. 



,+0.6 



7 kpc, and the lens star magnitude is 
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Fig. 1.12. The OGLE-2005-BLG-71 light curve showing the planetary anomaly near 
thepeak.The triple peak (two large symmetric peaks surrounding a small peak) indicates 
that the source passed three cusps of a caustic, the middle one being weak (insets). The 
interval between peaks (and so cusps) is At — 3 days, implying that the companion mass 
must be small. 



The light curve of the second planet discovered by microlensing, OGLE-2005- 
BLG-71Lb, is shown in Fig. 1.12 (Udalski et al. 2005). This was a moderately 



high magnification event that would have reached a maximum magnification of 
^4max ~ 42 if the lens star had no planets. Because of the d 1/d ambiguity 
discussed in § |1.4| this event has two models that explain the major features of the 
light curve quite well. Fig. |1.3| shows the magnification patterns for these models, 
and for the trajectory of the lens, which is nearly perpendicular to the lens axis, the 
light curves for these different models are quite similar. From Udalski et al. (20051 



the physically interesting parameters of the best fit models are t E = 
q = 7.1±0.3x 10~ 3 , and d = 1.294±0.002 for the "wide" model and t E 



70.9 ± 3.3, 
= 73.9±3.5, 
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q = 6.7 ± 0.3 x 10~ 3 , and d = 0.758 ± 0.002 for the "close" model. However, the X 2 
difference between these two models is Ax 2 = Xylose — ^widc = 22.0, so the "wide" 
model is strongly preferred. 

OGLE-2005-BLG-71Lb was the first planet discovery with significant contri- 
butions from amateur astronomers, with critical observations near the two strong 
cusp approach peaks by Grant Christie of the Auckland Observatory and Jennie 
McCormick of the Farm Cove Observatory. 




-20 20 

days since 31.0 July 2005 UT 

Fig. 1.13. The observed light curve of the OGLE-2005-BLG-390 microlensing event and 
best-fit model plotted as a function of time. The data set consists of 650 data points from 
PLANET Danish (red points), PLANET Perth (blue), PLANET Canopus (Hobart, cyan), 
RoboNet Faulkes North (green), OGLE (black), and MOA (brown). The top left inset 
shows the OGLE light curve extending over the previous 4 years, whereas the top right 
one shows a zoom of the planetary deviation, covering a time interval of 1.5 days. The 
solid curve is the best binary lens model described in the text with q — 7.6 ± 0.7 x 10~ 5 , 
and a projected separation of d — 1.610 ± 0.008-Re. The dashed grey curve is the best 
binary source model that is rejected by the data, and the dashed orange line is the best 
single lens model. 

With a mass ratio of q = 7.1±0.3 x 10~ 3 , OGLE-2005-BLG-71Lb must certainly 
be a gas giant planet, but without further information such as measurement of finite 
source effects, the detection of the lens star or a measurement of the microlensing 
parallax effect, we cannot determine the properties of the host star or the planetary 
mass with much precision. Fortunately, we are able to detect lens star in a set of 
HST images, and the light curve yields weak detections of both a finite source size 
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and the microlensing parallax effect. So, we expect to determine the host star and 
planet masses and to convert their separation into physical units, but this analysis 
is not yet complete (Dong et al 2007, in preparation). 

The first low-mass planet discovered by microlensing was OGLE-2005-BLG- 
390Lb (iBeaulieu et al.l 120061), led by the PLANET Collaboration. This planet is 



currently tied with Gl 581c (Udry ej^al. 20071 as the lowest mass exoplanet orbiting 



a normal star yet to be discovered^ This event was detected through a planetary 
caustic deviation, and the amplitude of the deviation was significantly reduced by 
the finite angular size of the clump giant source star. If the planet were smaller by 
a factor of ~ 2, it would not have been detected in this event. As originally pointed 
out in Bennett & Rhie (19961 and discussed in § 1.3.3 a microlensing search for 



Earth-mass planets should focus on events with main sequence source stars. 




Fig. 1.14. Comparison of the OGLE-2005-BLG-390Lb planetary caustic (the black dia- 
mond shaped curve) with the source star size (red circle). The black line with the arrow 
show the motion of the source star. 



The OGLE-2005-BLG-390Lb light curve deviation does not show the charac- 
teristic features of a fold caustic crossing, like OGLE-2003-BLG-235, or of a cusp 
approach, like OGLE-2005-BLG-71. This is because the planetary caustic is smaller 
than the source star's angular radius of 8* = 5.3 ± 0.7/xas is larger than the plan- 
etary caustic, as shown in Fig. 1.14 Because the light curve does not show these 



characteristic binary-microlensing features, we must consider a non-planetary ex- 
planation for the light curve involving the lcnsing of a binary source star by a single 



star lens. Gaudi (19981. However, as the Fig. 1.13 shows, a binary source model is 



a poor fit to the data, as it fails to account for the Perth and Danish data near 
the end of the perturbation. Formally, the binary source model increases the fit \ 2 
by Ay 2 — 46.25 with one fewer degree of freedom. These data are also sufficient to 



The minimum mass of M p > 5.O3M0 is often quoted for Gl 581c, but the M p sini 
ambiguity of the radial velocity method implies that the median predicted mass is 
M p = 5.5M®. This is the appropriate number to compare to other detection methods. 
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avoid a possible degeneracy in the planetary parameters for such events that occurs 



when the wings of the deviation are poorly sampled (Gaudi & Gould 19971 
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Fig. 1.15. Bayesian probability densities for the properties of the planet and its host star: 
(a), the masses of the lens star and its planet (M, and M p respectively), (b), their distance 
from the observer Dl, (c), the three dimensional separation or semi-major axis a of an 
assumed circular planetary orbit; and (d), the orbital period Q of the planet. The bold, 
curved line in each panel is the cumulative distribution, with the percentiles listed on the 
right. The dashed vertical lines indicate the medians, and the shading indicates the central 
68.3% confidence intervals, while dots and arrows on the abscissa mark the expectation 
value and standard deviation. The medians of these distributions yield a M v = 5.5^ 2 ' 7 M^ 
planetary companion at a separation of d = 2.6 g AU from a Af* = 0.22 ^q'jJ Mq Galac- 
tic Bulge M-dwarf at a distance of Dl = 6.6±1.0kpc from the Sun. The median planetary 
period is Q — 9^? years. 



The microlensing model for this event directly determines the planetstar mass 
ratio, q = 7.6 ± 0.7 x 10~ 5 , the projected planetstar separation, d = 1.610 ± 0.008, 
the Einstein radius crossing time, t& = 11.03 ± 0.11 days, and the source radius 
crossing time, i* = 0.282 x 0.010 days. With the value for 9* mentioned above, this 
yields the angular Einstein radius, Be — 0.21 ± 0.03 mas, from eq. 1.16 and the 



mass-distance relation from eq. 1.17 This mass-distance relation can be combined 
with a standard Galactic model in a Bayesian analysis to estimate the probability 



distribution of the lens system parameters ( Alcock et al. 1995 1996 Poindexter et 



al. 2005 Dominik 20061. The results of such an analysis are shown in Fig. 1.15 



following the method of Dominik (20061, and nearly identical results are obtained 
using the Galactic model and mass function parameters of Bennett & Rhie ( 2002 1 . 
This analysis gives a 95% probability that the planetary host star is a main-sequence 
star, a 4% probability that it is a white dwarf, and a probability of,l% that it is a 
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neutron star or black hole. The median parameters shown in Fig. |1.15| imply that 
the planet receives radiation from its host star that is only 0.1% of the radiation 
that the Earth receives from the Sun, so the probable surface temperature of the 
planet is 50 K, similar to the temperature of Neptune. 




-2 




-0.4 -0.2 0.2 0.4 

At [HJD - 2453491.875] 



Fig. 1.16. Top: Data and best-fit model for OGLE-2005-BLG-169. Bottom: Difference 
between this model and a single- lens model with the same single lens parameters (to, uo, tE, 
and p). ft displays the classical form of a caustic entrance/exit that is often seen in binary 
microlensing events, where the amplitudes and timescales are several orders of magnitude 
larger than seen here. MDM data trace the characteristic slope change at the caustic exit 
(At — 0.092) very well, while the entrance is tracked by a single point at At = —0.1427. 
The dashed line indicates the time to- Inset: Source path through the caustic geometry. 
The source size, p, is indicated. 



As discussed in 
star is detected. However 



1.4.1 the lens star mass can be determined directly if the lens 
this will be quite difficult for OGLE-2005-BLG-390L, 



because the source is a giant star. For the median mass and distance to the lens 
system, the lens star would be fainter than the source by a factor of 2000 in the 
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-ftf-band. So, the detection of the lens star may require many years for the relative 
proper motion of // rc i = 6.8mas/yr and the development of new instruments for 
large ground-based or space telescope. 

OGLE-2005-BLG-169 was the third published event from the 2005 season and 



the second low-mass planet found by microlensing (Gould et al. 20061. This was 



a very high magnification event, with a peak magnification of A max ~ 800, and 
it's light curve is shown in Fig. |1.16| The bottom panel of Fig. |1.16| indicates that 
the planetary deviation has a maximum amplitude of about 4% compared to the 
light curve of the same event without a planet. Such low amplitude deviations 
characteristic of the very weak caustics due to low-mass planets near the Einstein 
ring. However, it is only part of the caustic curve that is so weak. If the source 
would have passed on the other side of the host star and crossed the backwards "C" 
shaped part of the caustic in the inset of Fig. |1.16[ the planetary signal would have 
been very much stronger. But, it order to detect the low amplitude signal due to the 
caustics actually crossed by the source star, it was quite helpful to have continuous 
observations over the course of three hours from the 2.4 m MDM telescope. 



star: M,/M G 




Fig. 1.17. OGLE-2005-BLG-169 lens property figure. Bayesian probability densities for 
the properties of the planet and its host star if it is a main sequence star, (a) The masses 
of the lens star and its planet (M, and M v respectively), (b) the separation, (c) their 
distance from the observer (Dl)\ and (d) the I-band brightness of the host star. The 
dashed vertical lines indicate the medians, and the shading indicates the central 68.3% 
and 95.4% confidence intervals. All estimates follow from a Bayesian analysis assuming 
a standard model for the disk and bulge population of the Milky Way, the stellar mass 
function of Bennett & Rhie (2002). 
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The analysis of Gould et al. (20061 indicated a super-Earth mass planet with 
M p = 13+gAf ffi orbiting a star of M* ~ 0.49M Q . Such a planet, like OGLE-2005- 
BLG-390Lb, would be invisible to other planet detection methods. High magnifica- 
tion events also place tight constraints on the presence of additional planets, and 
in the case of OGLE-2005-BLG-169, Jupiter-mass planets can be exclude from the 
separation range 0.6-18 AU and Saturn-mass planets can be excluded from the range 
1-11 AU. 

The precise masses of the host star and planet have not yet been determined 
because the host star has not been detected. Thus, the lens system properties can 
only be determined by a Baycsian analysis, as was done for OGLE-2005-BLG-390Lb 
in Fig. 1.15 This analysis uses the parameters from the microlensing light curve, 
including the Einstein radius crossing time of ts — 43 ± 1.3 days, the source radius 
crossing time of i* = 0.019 ± 0.01 days, and the lens-source relative proper motion 

Of /i re i 



8.4 ± 0.6mas/yr. The results of this analysis are presented in Figure 1.17 



These have assumed a Han-Gould model for the Galactic bar (Han & Gould 19951, 
a double-exponential disk with a scale height of 325 pc, and a scale length of 3.5 kpc, 

( [20021 . 
(20061, 



as well as other Galactic model parameters as described in Bennett & Rhic 



Because this model is different from the Galactic model used by Gould et al 



the resulting parameters differ slightly from their results. We find a lens system 
distance of Dj, = 2.7^g'g kpc, a three dimensional star-planet separation of a = 



3-3 ^g'g AU and main sequence stellar and planetary masses of M* 

1+5 



0.52 



+0.19 
-0.22 



Mr. 



and M p = 14 «M 9 . If we assume that white dwarfs have an a priori probability to 
host planets that is equal to that of main sequence stars (at the separations probed 
by microlensing), then there is a 35% probability that the host star is a white dwarf. 
The possibility of a brown dwarf host star is excluded by the light curve limits on 



the microlensing parallax effect ( Gould et al. 2006 1 



Fig. 1.17 'd) shows the probability distribution of the /-band magnitude of the 
planetary host star compared to the source star at I = 20.58 ± 0.10. The implied 
planetary host star brightness distribution has a median and 1-cr range of Ii cns = 



21.9 



+0.7 



but the most interesting feature of this figure is that the probability of a 



main sequence lens fainter than / = 23 vanishes. This is because the mass-distance 
relation, cq. |1.17| ensures that the lens star will be nearby and at least at bright as 
I = 23, even if it is at the bottom of the main sequence at Af* = 0.08Af Q . In fact, 
the microlensing parallax constraint from the light curve yields a lower limit for 
the lens star mass of M* ^ 0.14M Q . Thus, the planetary host star must be at least 
16% of the brightness of the combined lens plus source star blended image, and this 
implies that it will be detectable if it is not a stellar remnant. Plus, the relatively 
rapid relative proper motion, fi Te i = 8.4 ± 0.6 mas/yr, of OGLE-2005-BLG-169L, 



implies that the lens-source separation is already detectable with HST (Bennett et 



al. 2007a I, as discussed in 



1.4.1 



One of the most interesting consequences of the discoveries of OGLE-2005-BLG- 
390Lb and 169Lb is that super-earth planets of ~ 5-15M0 are likely to be quite 
common. Gould et al. (20061 combined these detections with null results from very 



high magnification events (Abe et al. 2004 Dong et al. 20061 plus samples of 
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lower magnification events (Albrow et al. 2001 Gaudi et al. 20021 to solve for the 
fraction, / se , of stars with planets of mass ratio ~ 8 x 10 -5 at the separations of 
1.5-4 AU, where microlensing is most sensitive. They found that the median and 
90% confidence level upper and lower limits are / se = 0.38 ^'22 > based on the two 
planets discovered and the accumulated null results. The 90% c.l. lower limit is 
fse > 16%. This is significantly higher than the fraction of F, G, and K stars with 
Jupiter-mass planets in this 1.5-4 AU separation range. This fraction of stars with 
Jupiters at this separation can be estimated from Butler et al. (2006) to be fj ~ 3%. 



Thus, these cool, super-earth planets appear to represent the most common type 
of exoplanet yet discovered. This would seem to confirm the prediction of the core- 
accretion theory that ~ 10M® planets form much more frequently than gas giants, 



like Jupiter (Ida & Lin 2004 Laughlin, Bodcnhcimcr & Adams 2004), although 



this may not be incompatible with the disk instability theory ( Boss 2006 1 
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Fig. 1.18. Two views of the OGLE-2006-BLG-109 light curve, which is the first multi- 
planet microlensing event with a planet of slightly less than a Jupiter mass (q = 1.35 x 10 -3 
at ~ 2.7 AU and a planet of slightly more than a Saturn mass (q = 4.9 x 10" 4 ) at ~ 5.4 AU. 
The signal is dominated by the Saturn-mass planet because it is close to the Einstein ring 
at d — 1.4, and there are two pairs of caustic crossing features (at t — 3822.5, 3822.9 and 
t = 3830.2, 3831.2) and a cusp approach (at t = 3834.1) due to the Saturn-mass planet. 
The Jupiter-mass planet planet is further from the Einstein ring at d = 0.63, so its signal is 
limited to the highest magnification part of the light curve and is responsible for the cusp 
approach feature at t — 3831.65. Both planetary orbital motion and microlensing parallax 
must be included to obtain an acceptable model for this event. 



The final event that we will present is OGLE-2006-BLG-109, which is much more 
complicated than the other events (Gaudi et al. 2007; Bennett et al. 2007, both in 
preparation). The light curve for this event is shown in Fig. 1.18 while the central 
caustic configuration is shown in Fig. |1.19| This is the first microlensing event with 
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Fig. 1.19. The central caustic OGLE-2006-BLG-109 configuration for is shown at 3-day 
intervals from t — 3820 (shortly before the first caustic crossing) through t — 3835 (a day 
after the final cusp approach) . The time-order of the different color caustic curves is red, 
magenta, green, black, cyan, blue. The grey curve is the source trajectory, which is curved 
due to the microlensing parallax effect (i.e. the orbital motion of the Earth) and the small 
circle that the source trajectory in the left, close-up panel shows the source star radius. 

two detected planets, and it also shows clear signals of planetary orbital motion and 
microlensing parallax. These effects are detectable because the Saturn-mass planet 
has a projected separation that is close to the Einstein ring, which causes its caustic 
to become quite extended, fts effects are visible for 11 days. 

Another notable feature of OGLE-2006-BLG-109 is that the lens is > 5 times 
brighter than the source. It is detectable (although not completely resolved) in the 
best seeing (0.7") OGLE images and is clearly visible in K and f/-band adaptive 
optics images from the Keck telescope. As a result, there are two methods to deter- 
mine the lens star mass: the combination of the 8e determination from the finite 
source effects and the microlensing parallax effect yields the lens mass via eq. |1.19[ 
while the lens star detection give the lens mass with the help of mass-luminosity 
relations, as discussed in § |1.4.1| However, one complication is that there is some 
degeneracy in the effects of microlensing parallax and the planetary orbital motion 
on the microlensing light curve. On the other hand, the planetary orbital motion 
parameters yield information about the orbits that haven't been detected before 
in a microlensing event. So, this event will yield much more information about the 
OGLE-2006-BLG-109L planetary system than was anticipated for any planetary 
microlensing event. 

1.6 Future Programs 

Our experience with the existing microlensing planet search programs provides in- 
dications of how the sensitivity of future microlensing surveys can be improved. At 
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present, the OGLE and MOA groups are each able to independently discover more 
than 500 microlensing events per year. There is a great overlap between the discov- 
eries of these two groups, but the total number of events discovered every year is 
probably about 700. This is at least an order of magnitude larger than the global 
follow-up groups can hope to follow. The follow-up groups optimize their observa- 
tions by focusing on high magnification events. However, many of the shorter time 
scale high magnification events are not recognized as such in time, and so a large 
fraction of the high magnification events are not searched for planets. 

The solution to this problem is to observe many microlensing events in each 
image with a global network of very wide FOV telescopes that can observe 10 
square degrees or more of the Galactic bulge at 15-20 intervals. The new 1.8m 



MOA-II telescope (Hearnshaw et al. 20051 with a 2.2 square degree FOV CCD 



camera (Yanagisawa et al. 20001 that began operation in 2006 is the first telescope 
that meets this requirement, and the OGLE group plans to upgrade to a 1.4 square 
degree OGLE-4 camera in time for the 2009 Galactic bulge observing season. With 
MOA-II in New Zealand, and OGLE- IV in Chile, all that is needed is a very wide- 
FOV microlensing survey telescope in Southern Africa. A number of groups are 
pursuing funding for such a telescope. 
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Fig. 1.20. The number of planet detections expected per year as a function of planet 
mass is shown for proposed future space and ground-based microlensing surveys under the 
assumption of one planet per star in the indicated separation ranges. The space-based sur- 
vey has its most significant advantage over the ground-based survey at separations smaller 
(0.5-1.5 AU) and larger (5-15 AU) than the Einstein radius, because a space-based survey 
is able to resolve bulge main sequence stars and detect moderate amplitude planetary 
signals when the magnification due to the stellar lens is small. 
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Simulations of such a system have been performed by Bennett ( 2004 1 and Gaudi 
(2007, private communication), and estimates of the sensitivity of the global network 
consisting of MOA-II, OGLE-IV and an OGLE-IV-like system in South Africa are 
presented in Fig. |1.9| and Fig. |1.20| The improvement in sensitivity with such a 
network in the mass vs. semi-major axis plane is shown in Fig. 1.9 with the light 
and dark red curves showing the sensitivity of the current surveys and the future 
very wide-FOV network, respectively. This network will extend the sensitivity of 
the microlensing method down to an Earth mass at planet-star separations close to 
the Einstein ring radius (~ 2-3 AU) . 

The separation range where ground-based microlensing is most sensitive, 1-5 AU 
corresponds to the vicinity of the so-called "snow-line" which is the region of the 
proto-planetary disk where it is cold enough for water-ice to condense, The density 
of solids in the proto-planetary increases by a factor of ~ 4 across the "snow- 
line" and as a result, the core accretion theory predicts that this is where the most 



massive planets will form (Ida & Lin 2004 Laughlin, Bodenheimer & Adams 2004 



Kennedy et al. 2006 1. According to this theory, giant planets form just outside the 
"snow line" where they can accrete ~ IOMq of rock and ice to form a core that 
grows into a gas giant like Jupiter or Saturn via the run-away accretion of Hydrogen 
and Helium onto this core. However, this theory also predicts that the Hydrogen 
and Helium gas can easily be removed from the proto-planetary disk during the 
millions of years that it takes to build the rock-ice core of a gas-giant. Thus, if the 
core accretion theory is correct, rock-ice planets of ~ 10M® that failed to grow into 
gas giants should be quite common, although it is possible to form such planets in 
the competing gravitational instability theory (Boss 20061. 

The number of planetary microlensing event detections expected per year is 
shown in Fig. 1.20 assuming an average of one such planet per star, with conservative 
assumptions regarding photometric precision. The assumption of an average of one 



such planet per star is certainly too optimistic for Jupiter mass planets (Gaudi et 



al. 


2002 


Butler et al. 


2006) 



2006), but it is closer to reality for super-earths, like OGLE- 



2005-BLG-390Lb and OGLE-2005-BLG-169Lb (Beaulieu et al. 2006 Gould et al. 



2006 I . It could very well be accurate for Earth-mass planets where the weaker two- 



body gravitational interactions allow two planets to orbit in the separation range 
corresponding to the bins in Fig. 1.20 (Our own solar system is an example of this.) 

Another future development that is already funded is a global network of robotic 
telescopes dedicated to monitoring transient events like planetary microlensing 



events, known as the Las Cumbres Global Telescope Network (Brown et al. 2007). 



Ideally, this network would routinely observe high magnification microlensing events 
and planetary deviations discovered in progress with an very high cadence, such as 
that provided by the MDM telescope for OGLE-2005-BLG-169 (see Fig. |1.16| This 
would enable the very wide-FOV survey telescopes to maintain their normal sam- 
pling strategy so that other planetary microlensing events would not be missed. 
This might add to the planet detection efficiency substantially, but such a system 
is more difficult to model. 
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1.6.1 The Ultimate Exoplanet Census: Space-Based Microlensing 



The ultimate census for virtually all types of exoplanets would be a space-based 
microlensing survey (Bennett & Rhie 2002 Bennett et al. 2007b I. Such a survey 



could provide a statistical census of exoplanets with masses > 0.1M® and orbital 
separations ranging from 0.5 AU to oo. This includes analogs to all the Solar Sys- 
tems planets except for Mercury, as well as most types of planets predicted by planet 
formation theories. This survey would determine the frequency of planets around 
all types of stars except those with short lifetimes. Close-in planets with separations 
< 0.5 AU arc invisible to a space-based microlensing survey, but these can be found 
by Kepler (Basri et al. 20051. Other methods, including ground-based microlens- 



ing, cannot approach the comprehensive statistics on the mass and semi-major axis 
distribution of extrasolar planets that a space-based microlensing survey will pro- 



vide. Detailed simulations of a space-based microlensing survey (Bennett & Rhie 



20021 have been used to determine the sensitivity of such a mission, and Figs. 1.9 



and 1.20 show the sensitivity of the proposed Microlensing Planet Finder (MPF) 



mission (IBennett et all [2004|. These fi gures also show that the sensitivity of a 



ground-based microlensing survey to terrestrial planets is limited to the vicinity of 
the Einstein radius at 2-3 AU. This is because ground-based survey generally re- 
quires moderately high magnification A ^ 10 in order to resolve the source star well 
enough to get the moderately precise photometry that is needed to detect planets 
with the microlensing method. A space-based microlensing survey would generally 
resolve the source stars, so planets further from the Einstein radius can be detected 
via their light curve perturbations at relatively low magnification from the lensing 
effect of the planetary host star. 




Percent Uncertainty Percent Uncertainty 



Fig. 1.21. (a) The simulated distribution of stellar masses for stars with detected ter- 
restrial planets. The red histogram indicates the subset of this distribution for which the 
masses can be determined to better than 20%. (b) The distribution of uncertainties in 
the projected star-planet separation, (c) The distribution of uncertainties in the star and 
planet masses. Note that it is the two-dimensional projected separation that is presented 
here, and we have not included the uncertainty in the separation along the line-of-sight as 
was done in Figure [T.17| 



1 Detection of Extrasolar Planets by Gravitational Microlensing 



37 



A space-based microlensing survey is also able to detect most of the planetary 
host stars for most planetary microlensing events. Using the methods described in 
and in more detail in Bennett et al. (2007a I, this allows the determination 



1.4.1 



of the star and planet masses and separation in physical units. This can be accom- 
plished with HST observations for a small number of planetary microlensing events 
(Bennett et al. 2006), but only a space-based survey can do this for hundreds or 
thousands of planetary microlensing events that future surveys would expect to 
discover. Fig. |1.17| shows the distribution of planetary host star masses and the 
predicted uncertainties in the masses and separation of the planets and their host 



stars (Bennett et al. 2007a I from simulations of the MPF mission. The host stars 



with masses determined to better than 20% are indicated by the red histogram in 
Fig. 1.17[ a), and these are primarily the host stars that can be detected in MPF 
images. Ground-based microlensing surveys also suffer significant losses in data cov- 
erage and quality due to poor weather and seeing. As a result, a significant fraction 
of the planetary deviations seen in a ground-based microlensing survey will have 
poorly constrained planet parameters due to poor light curve coverage (Peale 2003 1. 
(These poorly characterized detections are not included in Fig. 1.20 however.) 

Proposed improvements to ground-based microlensing surveys can detect Earth- 
mass planets in the vicinity of the "snow-line" which is critical for the understanding 
of planet formation theories (Gould et al. 2007). But such a survey would have its 
sensitivity to Earth-like planets limited to a narrow range of semi-major axes, so 
it would not provide the complete picture of the frequency of exoplanets down to 
0.1M® that a space-based microlensing survey would provide. Such a survey would 
probably not detect the planetary host stars for most of the events, and so it cannot 
provide the systematic data on the variation of exoplanet properties as a function 
of host star type that a space-based survey will provide. 

A space-based microlensing survey, such as MPF, will provide a census of ex- 
trasolar planets that is complete (in a statistical sense) down to O.IM^ at orbital 
separations > 0.5 AU, and when combined with the results of the Kepler mission 
a space-based microlensing survey will give a comprehensive picture of all types of 
extrasolar planets with masses down to well below an Earth mass. This complete 
coverage of planets at all separations can be used to calibrate the poorly understood 
theory of planetary migration. This fundamental exoplanet census data is needed 
to gain a comprehensive understanding of processes of planet formation and migra- 
tion, and this understanding of planet formation is an important ingredient for the 
understanding of the requirements for habitable planets and the development of life 



on extrasolar planets (Bennett et al. 2007b I 



The basic requirements for a space-based microlensing survey are a 1-m class 
wide field-of-view space telescope that can image the central Galactic bulge in the 
near-IR or optical almost continuously for periods of at least several months at a 
time. This can be accomplished as a NASA Discovery mission, as the example of 
the MPF mission shows, but it could also be combined with other programs that 
require an IR-optimized wide-FOV space telescope, as long as a large fraction of 
the observing time is devoted to Galactic bulge observations. As Fig. |1.9| shows, 
there is no other planned mission that can duplicate the science return of a space- 
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based microlensing survey, and our knowledge of exoplanets and their formation 
will remain incomplete until such a mission is flown. 

Thus, a space-based microlensing survey is likely to be the only way to gain a 
comprehensive understanding of the nature of planetary systems, which is needed 
to understand planet formation and habitability. The proposed Microlensing Planet 
Finder (MPF) mission is an example of a space-based microlensing survey that can 
accomplish these objectives with proven technology and a cost that fits comfortably 
under the NASA Discovery Program cost cap. 
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